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Abstract 

We show that the only Tolman models which permit a Vaidya limit 
are those having a dust distribution that is hollow — such as the 
self-similar case. Thus the naked shell-focussing singularities found 
in Tolman models that are dense through the origin have no Vaidya 
equivalent. This also casts light on the nature of the Vaidya metric. 
We point out a hidden assumption in Lemos' demonstration that 
the Vaidya metric is a null limit of the Tolman metric, and in 
generalising his result, we find that a different transformation of 
coordinates is required. 
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Introduction 

Recently Lemos [1992] showed that the Vaidya metric, describing radially directed inco- 
herent radiation (spherically symmetric null dust), can be obtained from the Tolman metric, 
which represents a spherically symmetric distribution of pressureless matter (dust), by taking 
a null limit. This surprising and intriguing insight was inspired by the very strong similarities, 
quantitative as well as qualitative, between the naked shell-focussing singularities (discovered 
by Eardley and Smarr [1979]) that appear in the self-similar forms of these metrics at the mo- 
ment the crunch singularity forms [see for example Christodoulou 1984, Newmann 1986, Ori 
and Piran 1987, Rajagopal and Lake 1987, Hellaby and Lake 1988, Waugh and Lake 1988, 1989, 
Grillo 1991, Lemos 1991, and extensive references in footnote 2 of Lake 1992]. We discuss the 
nature of an origin of spherical coordinates in Tolman models, and show that a Vaidya limit 
cannot be extended to such a point. We show that one of Lemos' assumptions can be relaxed 
if a different coordinate transformation is used. 

The incoming Vaidya metric [Vaidya 1951, 1953, see also Lindquist Schwarz and Misner 
1965] is 

ds 2 = 2dvdR- (i-^r)^ 2 + R 2 dSl 2 (1) 

where dVt 2 = d6 2 + sin 2 (0)d0 2 is the metric on a 2-sphere, the areal radius is positive, R > 0, 
and M = M(v) > is an arbitrary function of the null coordinate v, representing the effective 
gravitational mass inside v. The only non-zero Einstein tensor component and the Kretschmann 
scalar K = R a ^ 5 R al3 ^s are 

G v = —M* (2) 

vv R2 

- ^ (3) 

where * = Jjj, and superscripts V and T are used where necessary to distinguish quantities in 
the Vaidya and Tolman models. 

The Tolman metric [Lemaitre 1933, Tolman 1934] uses synchronous coordinates that are 
comoving with the dust particles, 

ds* = _ df 2 + -^jdr 2 + R 2 dil 2 (4) 

where ' = J^, / = /(r) is an arbitrary function of coordinate radius that determines the local 
spatial geometry, as a function of r [see Hellaby and Lake 1985, Hellaby 1987]. The areal radius 
R = R(t, r) is a solution of 

B^'f.f (5, 

where ' e |, and M = M(r) > is a second arbitrary function. Comparing this equation 
with its Newtonian analogue for the kinetic plus potential energy of a radially moving particle 
of mass m at a distance x from the centre of a spherically symmetric dust cloud with density 
distribution pn(x) 



where 



X 



M N (x) = / 4irx 2 p N (x)dx (7) 
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we obtain the interpretation that M(r) is the gravitational mass within co-moving radius r, 
and f(r) is twice the energy per unit mass of the particles at r. (The principal difference 
between these two equations is the replacement of the radial distance x by the areal radius R.) 
For / > (or rather Rf/M > 0) the evolution of the areal radius for a collapsing model is 
hyperbolic 



R 

(sinh?7 — if) 



M 

T 



(cosh T] 



M 



(8) 
(9) 



where the third arbitrary function a = a(r) gives the time at which R = — the big crunch. 
(Parabolic and elliptic solutions exist for / = and / < 0.) Since the pressure is zero, the dust 
particles (which stay at constant r, 9, 0) follow geodesies of the spacetime. It can be shown in 
general [Hellaby and Lake 1984, 1985] that for the collapsing models 



-SB)- 



The density and the Kretschmann scalar are given by 

2M' 




R 



K 1 



R 2 R' 
48M 2 



32MM' Y2M 



R 5 R' 



R*R' 2 



(10) 

(11) 
(12) 



Lemos' Method 

We here outline the approach used by Lemos, although we find it convenient to delay taking 
the null limit until a slightly later stage in the working. He initially makes the assumption of 
self-similarity in both metrics, for simplicity, and later drops it. That assumption is not made 
here. The Tolman line element @ may be transformed from coordinates (t, r) to (t, R) by 
means of 

dR = Rdt + R'dr -> R'dr = dR- Rdt (13) 

which leads to 

l2 / 2M\ dt 2 2R , Iri dR 2 n2 irV2 

where the new g tt has been simplified using 

Now the limit of interest is that in which / is allowed to diverge, while M and R are both 
required to remain finite 

/^+oo, < R,M < oo (15) 

Eq (H) shows that in this limit coshrj must also diverge, so that coshr/ — > smlir? — > e v /2 and 
(H) plus (H) simplify to 

R^yff(a-t) (16) 
and, for finite R, (a — t) must be vanishingly small. Similarly ([5]) for collapsing models becomes 



fl--V/ ( 17 ) 
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and the derivative of fll6|) (or alternatively substituting for (a — t) and R from (|16|) and fll?!) in 
(H)) gives 

Rf . r~ 



t>i R f , i r f 

R —r-r + a v J 



2/ 



Lemos then states that the transformation 

t 



v 



R 

7f + l 



(19) 



converts (fH| ) into the Vaidya metric (|1|), in the limit / — > oo. Since (|19|) and (|T6| ) imply that 



(20) 



the new coordinate becomes a function of r only, in the limit, so we can now write M — > M(u). 

We note however that a constant /, inherited from the self-similar case, must still be 
assumed in order to get this result. If we don't make this assumption, then ( |19"D leads to 



dt dR ( t R 

" 77 + T~\2F~ 2 + T 2 , 

dt dR 



f'dr 



t R\ ,{dR-Rdt) 



VT f V2/ 3/2 PJ 

and, after substituting for t, R, and R' from (0)-(0), to 



./ 



R' 



(21) 
(22) 



dv 



dt dR" 

77 + T, 



[l-X) 



where 



X 



R + a^/J 



R + 2p/\a'/f) 



so that (O) in the limit becomes 



ds 2 



f 



1 + /, 



/ 



dR 2 + 



f 



i + fj Vi-x 



i +7 u 



2M 



./ 



i + vi-x 



i 



i - 



2M 
~R 



dv 1 + i? 2 rffi 2 



(23) 
(24) 

dvdR 
(25) 



The limiting behaviour of X is not at all clear, as the relationship between a(r) and /(r) is 
arbitrary in general, and the limiting behaviour of a is not specified. 



The Problem of the Origin and the Form of f(r) 

In all Tolman models describing a collapsing dust cloud which exhibit a naked singularity, 
this singularity occurs at the moment of collapse £ = a, at the origin (r = being the natural 
choice). The origin of spherical coordinates is specified by R(t, r = 0) = 0, V t and we also have 
R(t,r — 0) = 0, V t, which, by eq (|), requires M(0) = as well as /(0) = at the origin, for 
example the homogeneous case (dust FLRW). Can we extend Lemos' result for the null limit 
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to cases where / does not diverge at the origin? Clearly the functional form of f(r) must allow 
/(0) = — for example / = pr 2 , p — > oo. 

Consider cases with / finite at r = 0, such as the non-parabolic self-similar models. As- 
suming R,M > 0, eq (J^j shows that R(t, r = 0) > / ^ 0. If r = is approached along 
a constant t surface, with a(r) finite near r = 0, (§) and (|9D show that either (a) M — > 
so that 1] — > oo and i? — > yfft, or (b) M remains finite so that rj and R also remain finite. 
Case (a) represents a hollow dust cloud — it can be matched at r = onto a vacuum Tolman 
(Minkowski) spacetime with M(r) = and a true origin at some negative r value where / = 0. 
Case (b) either (i) contains more dust inside r = 0, with the true origin again at / = — i.e. 
r = is not the centre of the cloud, or (ii) it contains the Schwarzschild vacuum inside r = 0, 
with no origin, or (iii) it contains a dust filled version of the Schwarzschild-Kruskal-Szekeres 
topology [Hellaby 1987]. In (ii) and (iii), / must pass through zero and reach —1 in order to 
form the throat, rising to / > in the second sheet, and M,R < oo everywhere that / < 0. 
Clearly particle worldlines having /(0) > are not at the origin, but they do collapse to zero 
and begin the formation of the singularity. 

Furthermore, since shell-focussing singularities do form in Tolman models with normal ori- 
gins, can the detailed similarity between the naked singularities of the two metrics be extended 
to such cases, or is constant / required? 

Note also that the coordinate r is eliminated by the first transformation ( |T5j ) and then 
effectively re-introduced, in the limit v = a(r)/y / /(r), via the second one flT9"|). Since the 
Tolman coordinate r is co-moving with the dust particles, and the Vaidya coordinate v is co- 
moving with the shells of radiation, one might expect v to be the direct limit of r. This is 
consistent with the interpretation of / as an energy parameter which goes to infinity, meaning 
that the Tolman particle geodesies are asymptotically null. Since a particle staying at the origin 
of spherical symmetry cannot be moving at light speed, this suggests that a Vaidya limit is not 
achieveable here. 



The Null Limit for General / and a 

Consider approaching the origin on a constant rj surface. Equation (^) shows that Rf /M 
remains constant, whereas (||) shows that both M/R and / go to zero there. Thus the Vaidya 
limit could be described by 

^-^oo, < R,M < oo (26) 



which doesn't necessarily require / — » oo at r = 0. The limiting forms (|I6| ) -([T8"| ) of R, R, and 
R' are unchanged by this adjustment. 

Starting again from (|J), we transform from (t, r) to (R,r) as our coordinates, thus substi- 
tuting for t rather than r, 

dR = Rdt + R'dr -> dt = (dR - R'dr) / R (27) 

and apply (W) to simplify the resulting g rr 



ds 2 = -± d R 2 + 2^-dRdr - ( 1 - —) ?—^dr 2 + R 2 dQ 2 (2* 

R 2 R 2 V R J (1 + f)R* V 



From ( |P7| ) and (|I8D , we have the following limiting forms of the extra factors that don't appear 
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in® 

J_ 1 B!_ (_a/_ Rf_) R' 2 ( f \j a' Rf 1 2 

' ^^l77 + 2pj ' (i + f)R 2 ~* vTT7y 177 + 2f j ( } 

The limiting transformation takes care of the first term in the curly brackets (a'/y/J), but 
not the second (Rf'/2f 2 ), and without knowing both a(r) and f(r) — i.e. f(a) — it can't be 
discounted. The second term is dominant if 

(Rf'\ ( a'\ R df 

W) ' \7lJ = W^Ta - °° (30) 

An example of an f(a) that makes the second term dominant almost everywhere is 

/ = a In (p) + sin {p n a) , n const. , p — > oo (31) 
df/da f + ae nf/a cos(ae n//a ) 



f3/2 a p/2 



(32) 



but this wildly oscillating form is very unrealistic. The conditions for no shell crossings [Hellaby 
and Lake 1985] for collapsing hyperbolic Tolman models require /' > and a' > wherever 
M' > 0, i.e. df/da > 0, so adding a linear term to remove negative gradients gives a vanishing 
second term 

/ = 2p n a + sin{p n a) , n const. , p — > oo (33) 
df/da 2 + cos(//2) 



fV 2 2av7 



(34) 



The most rapid uniform divergence of df/da relative to / we have been able to find for df /da > 
still leaves (df / da) / f 3 ^ 2 vanishing. It is expressed in terms of computer notation " for raising 
to the power, 

/ = ap"(p~(p\..(p"a))) , p^co (35) 
df/da ln(/)lnln(/)lnlnln(/)... 



PI 2 ay/7 



(36) 



However, at a single point (or a finite number of discrete points) the divergence behaviour can 
always be made arbitrarily rapid, e.g. 



/ = pa + a pn , n const. , p — » oo (37) 
df/da 1 

-> f n - 3/2 , a = 1 (39) 



< a < 1 (38) 



Hf) 
In(o)v7 



a > 1 (40) 



Consequently, we now introduce the following transformation, 



f r a' R 

" = L 7J dr ~v (41) 

( a' Rf'\ , <ifi 

* = {77 + w) ir -v (42) 
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which incorporates both terms in the brackets of (p9|), and which converts fl28|) to 



ds 2 



1 



4/(1 + /) 



2M\ 



f 



1 + /, 



1 



dR 2 4 
2M 



R 



(1 + /) 
cfo 2 + #W 



2M\ 
~RJ 



dvdR 



(43) 



It is already clear from (p29|) as well as this equation that / — > oo is indeed required everywhere 
to obtain the Vaidya metric as the limit. 



ds 2 — > 2dvdR 
(The alternative transformation 

=c?r — 



1 



T+7 



o v/r+7 



/ 



R 



dv 



2M\ 
~RJ 



/ 



dv 2 + R 2 dn 2 



(44) 



l + ZVv 7 / 2/2 



dr — 



'1 + / 



/ 



dR 



(45) 



does not succeed in removing the factor of //(l + /), and also leads to the wrong limit.) In the 
limit (15) then, it is evident from (|4i~|) that v becomes a function of r only, so that M — > M(v) 
holds once again. No assumptions about the functional form of / or the limiting behaviour of 
a were made to obtain the Vaidya metric as the null limit, and we find that the second term 
of pT| ) becomes negligible, even if the second term in the brackets of ( B2p doesn't. The new 
transformation (41) can also be re-written in the limit as 

dv 



dt dR 

77 + V 



(46) 



in order to recover ([!]) from fll~4|). Equations ([46]) and (f41"D are the revised versions of (^3|) and 



The overall transformation from Tolman to asymptotically Vaidya coordinates then is 



r a '( r Xdr- R{t ' r) 



f(r) 



2/(r) 



R 



R(t,r) 



(47) 
(48) 



where R(t,r) is given by Using the following limiting values of two of the partial 

derivatives of the inverse transformation 



dr 
dv 



f_ 
R> 



dt 
dv 



(49) 



R 



the Kretschmann scalar and the density may be converted to their Vaidya forms. Thus 



M* 



dM 




dM 


dr 


dv 


R 


dr 


t dv 



f 

M — 
R> 



i.e. 



R> f 



(50) 



ensures that the last two terms on the right of (1T2|) vanish, giving (|3|) in the limit. For the 
'density', (Q) is obtained from (ITTJ) in the limit by writing 



G 



v 



m 

dv 



f 



2M' 2 
R 2 R' = & 



M* 



(51) 
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The strengths of singularities are variously defined by [e.g. Tipler Clarke and Ellis 1980, Clarke 
and Krolak 1986] 

^ G = lim x ^\ 2 G a pk a k p or * R = lim x ^ Q \ 2 R a(3 k a k p (52) 

where k a is the tangent vector to a null geodesic with parameter A that hits the singularity at 
A = 0. From the above, and since ^ is a scalar, it is clear that the strengths of the Vaidya 
singularity, as measured along radial geodesies are given by the limits of the corresponding 
Tolman expressions. 



Conclusions 

Lemos originally demonstrated that the Vaidya model is a null limit of the Tolman model, 
by taking the limit / — > oo and assuming / = constant in this limit. His transformation 
was completely valid for models with constant /. However Tolman shell-focussing singularities 
also occur in models with matter at the origin. The existence of a normal origin of spherical 
coordinates at r = 0, (a — t) > in the Tolman model requires f(r = 0) = 0, and we have 
found this cannot be made consistent with a null limit. The Vaidya limit does indeed require 
/ — > oo, so it cannot be extended to a spherical origin, where /(0) = 0, or a Schwarzschild- 
Kruskal-Szekeres type topology, which requires / = — 1 in the throat. 

Thus we conclude that every Vaidya model is the limit of a hollow Tolman model, acquiring 
its arbitrary M(v) from a combination of M(r) and a(r), and must itself be hollow. If M(r = 
0) = 0, M(v = 0) = 0, then r = 0, v = is a collapsing shell of finite size surrounding 
Minkowski vacuum, and the limiting Vaidya model can form a shell-focussing. If Mir = 0) > 0, 
M(v = 0) > then it surrounds Schwarzschild vacuum, and no shell-focussing forms. In this 
latter case, the shells of incoming radiation (having / divergent) cannot pass through the throat 
(where / = — 1) and must hit the future singularity first. A dust filled interior is not possible in 
the limit, since a coordinate line cannot be co- moving with both a dust particle and a light ray, 
but it may be possible to have an intervening vacuum region. Since t = a on the singularity, 
and R is only finite on a collapsing shell of radiation where (a — t) is infinitesimal, the radiation 
is all at infinite R for any finite value of (a — t). 

If we assume that M* is finite, then it is apparent from (|41|) that a finite change in M 
and v requires an infinite change in a. It is interesting to note that a collapsing, unbound 
(i.e. hyperbolic) dust cloud of finite total mass may also have /, a — > oo and M finite in the 
asymptotic regions. At finite (a — t), i? is infinite, from (U)-(§), but as these particles collapse 
towards the crunch, R becomes finite when (a — t) is infinitesimal, and the Vaidya limit is 
achieved. In terms of Tolman time, this is infinitely long after the initial formation of the 
singularity, but only a finite retarded time in Vaidya coordinates. 

The new coordinate transformation (fi7])-(fi8|) — or fl27|) and fl42f) — makes no assumptions 
about the 3 arbitrary Tolman functions /, M, and a, and in particular the relationship between 
/ and a, beyond those normally made for a general Tolman model, and the limit / — > oo. 
Several important physical quantities — the Einstein tensor, the Kretschmann scalar, the null 
geodesies, and the strenghts of the singularities — all have the correct limit. This generalises 
Lemos' unification of the two metrics and their shell-focussing singularities. 
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